The Observability of Metal Lines Associated with the Lyman-alpha Forest 
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ABSTRACT 

We develop a prescription for characterizing the strengths of metal lines associated 



> 

O . 

, with Lya forest absorbers (LYFAs) of a given neutral hydrogen column density A^hi and 

metallicity [0/H]. This Line Observability Index (LOX) is line-specific and translates. 



I for weak lines, into a measure of the equivalent width. It can be evaluated quickly 



for thousands of transitions within the framework of a given model of the Lya forest, 
providing a ranking of the absorption lines in terms of their strengths and enabling 



Q^' model builders to select the lines that deserve more detailed consideration, i.e. those 

Q . that should be detectable in observed spectra of a given resolution and signal-to-noise 

^ , ratio. We compute the LOX for a large number of elements and transitions in two 

I cosmological models of the Lya forest at z ~ 3 derived from hydrodynamic simulations 

of structure formation. We present results for a cold dark matter universe with a 
^ cosmological constant; an = 1 cold dark matter model yields nearly identical results, 

' and we argue more generally that the LOX predictions are insensitive to the specific 

^ choice of cosmology. We also discuss how the LOX depends on redshift and on model 

parameters such as the mean baryonic density and radiation field. 

We find that the OVI (1032 A, 1038 A) doublet is the best probe of the metahicity 
in low column density LYFAs (Ahi ~ lO^^'^cm"^). Metallicities down to [0/H] ~ 
-3 yield OVI absorption features that should be detectable in current high-quality 
spectra, provided that the expected position of the OVI feature is not contaminated 
by HI absorption. The strongest transitions in lower ionization states of oxygen are 
OV(630 A), OIV(788 A), and 0111(833 A). These absorption hues are ah predicted 
to be stronger than the OVI feature, but even at redshifts 3 — 4 they will have to be 
observed in the ultraviolet, and they are extremely difficult to detect with present 
UV instruments, such as the Space Telescope Imaging Spectrograph (STIS). At lower 
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redshifts, detection of these lines may be possible in STIS spectra of the very brightest 
QSOs, while one may have to wait for next-generation instruments such as the Cosmic 
Origins Spectrograph (COS) to detect such lines in a number of high-redshift QSOs. 

The strongest metal lines with restframe wavelength larger than 912A associated 
with higher column density LYFAs at z ~ 3 are CIII (977 A) and Silll (1206.5 A), 
which peak at A^^hi ~ lO-'^^cm"^. Of the lines with rest wavelengths > 1216A, which 
can potentially be observed redwards of the Lya forest, the CIV(1548,1551) doublet 
is expected to dominate in all LYFAs, regardless of the value of A^hi- We argue that 
CIV and CII absorption may peak in different spatial regions, and that comparison of 
single-phase models of the CIV/CII ionization ratios with observed CIV/CII column 
density ratios can lead to an overestimate of the ionization parameter in the central 
parts of the absorbers. 

Subject headings: Intergalactic medium: ionization, quasars: absorption lines 

1. Introduction 

Over the past decade, observations have demonstrated that strong Lya forest absorbers 
(LYFAs) generally show associated metal line absorption (Meyer & York 1987; Lu 1991; Cowie et 
al. 1995; Womble et al. 1995; Songaila & Cowie 1996, hereafter SC). SC find CIV absorption in 
nearly all LYFAs at z ~ 3 with neutral column densities A^hi ~ 1-6 x 10"*^^ cm^^ and in ~ 75% 
of systems with A^hi ^ 3.2 x 10^^. The observed CIV column densities are consistent with the 
absorbers having a mean metallicity [C/H] ~ —2.5 and an intrinsic scatter in metallicity of about 
an order of magnitude (Ranch et al. 1996; Hellsten et al. 1997). An important question is whether 
or not a chemically pristine or extremely metal poor population of LYFAs exists, and if so, what 
the characteristic HI column densities of this population are. For logAni ^ 14.5 the associated 
CIV lines are close to the detection limits of SC, and it cannot yet be determined if these absorbers 
have the same metallicity distribution as those with higher A^hi- Some theoretical models of early 
metal enrichment predict that the mean metallicity declines steadily with decreasing A^hi (Gnedin 
& Ostriker 1997). 

Considerable progress has been made in the theoretical understanding of LYFAs. From the 
work of several groups (e.g. Cen et al. 1994; Zhang et al. 1995; Petitjean et al. 1995; Hernquist et 
al. 1996; Bi &; Davidsen 1997) a picture has emerged in which the LYFAs are interpreted in terms 
of "Gunn-Peterson" absorption in an inhomogenous IGM pervaded by a background ionizing 
radiation field, presumed to originate from QSOs and perhaps young, star-forming galaxies. 
This picture links the Lya forest directly to cosmological structure formation. The analysis 
of cosmological models of the LYFAs consists of producing artificial spectra by evaluating the 
absorption properties of the baryonic IGM along lines of sight through N-body and hydrodynamical 
simulations of structure formation. These spectra can then be analyzed in much the same way as 
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observed QSO spectra, and the resulting distribution functions in Nfji and linewidths are found to 
agree quantitatively with those observed, lending considerable support to the cosmological LYFA 
picture (e.g. Miralda-Escude et al. 1996; Dave et al. 1997; Zhang et al. 1997). 

The interpretation of the metal line data within the framework of these models provides a 
strong test of this picture. If a metal enrichment pattern of the baryonic IGM is specified, the 
properties of selected absorption lines are readily calculated from the knowledge of densities, 
temperatures, and UV radiation field along the lines of sight (Haehnelt et al. 1996; Ranch et 
al. 1996; Hellsten et al. 1997). Until now, the metal lines considered have been limited to 
OVI(1032,1038), NV(1239,1243), CIV(1548,1551), SiIV(1394,1403), and 011(1335). These lines 
have been chosen because they have already been identified with LYFAs in real QSO spectra. 
With the exception of the OVI doublet they all have rest wavelengths Aq > 1216 A, which means 
that for LYFAs at sufficiently high redshifts, they appear redward of the Lya forest and hence are 
relatively easy to detect. Instead of making an a priori selection of a few metal lines to include 
in a cosmological model for LYFAs, a more satisfactory approach would be to make the models 
predict which metal lines, out of thousands of candidates, should be observable (in spectra of a 
given resolution and S/N) in LYFAs of different column densities. Such an approach would allow 
a comprehensive and swift selection of lines that deserve a more detailed treatment, depending 
on the specific purpose of the modeling. It would provide a sharper test of the model and place 
more stringent constraints on the metallicities and UV radiation spectrum, by making it possible 
to compare the complete list of metal lines predicted to be observable to the list of lines actually 
detected in QSO spectra, as well as the observed relative strengths of these lines. 

This paper presents and applies a technique for implementing this comprehensive, a priori 
approach to metal line absorption predictions in cosmological models. In Section 2 we define a 
line observability index, which can easily be evaluated for a database of metal lines, within a 
particular model for LYFAs. This index is a measure of column density and hence allows the lines 
to be sorted by strength. It also predicts which lines should be detectable in LYFAs for a wide 
range of column densities and metallicities. In Section 3 we apply this technique to two specific 
cosmological models of the Lya forest, based on hydrodynamic simulations of structure formation. 
In Section 4 we discuss the dependence of the LOX on model parameters such as the mean baryon 
density and radiation field and estimate how it changes with redshift, and we argue that the 
LOX has only a very weak dependence on the underlying cosmology. We discuss the results and 
summarize the conclusions in Section 5. 



2. The Line Observability Index 

Detection limits for absorption lines can be conveniently expressed in terms of the equivalent 
width 
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where Fx/F^x denotes the received energy flux per wavelength interval, normalized to the 
continuum level. For absorption lines that are weak enough to be optically thin, the rest frame 
equivalent width Wxr is proportional to the column density of absorbing atoms along the line 
of sight 

2 

Wxr = ^NXr^f, (2) 

where / and A,, are the oscillator strength and rest frame wavelength of the line considered. The 
other symbols have their usual meaning. Lines that satisfy equation (|2[) are said to be on the 
linear part of their curve-of-growth. (For a derivation of this relation see, e.g., Spitzer 1978.) We 
will label an absorption line produced by an element Z in ionization stage i with rest transition 
wavelength A and oscillator strength / as Zlj- 

For a metal absorption line of metallicity [Z/H] = log (nz/nii) — log(n^/nH)Q associated 
with a LYFA of hydrogen column density A'hi we now define the following quantity: 

LOXiZij,[Z/H],Nni) = -17.05 + logiVni + 

which we will call the line observability index. The inspiration for this definition comes from 
rewriting equation (^) with the column density of the absorbing metal ions expressed in terms of 
A^Hi; [Z/H], and the solar ratio of the number densities of metal and H atoms {nz/nii)Q. The 
model-dependent ionization corrections are contained in the last term, which is a mean value 
of the ratio of the ionization fractions = n^/n^ of and in the region where the Lya 
absorption arises. Equation (^) assumes units of cm~^ for A'^hi and A for A; the choice of additive 
constant then implies LOX = log {Wrx/ 1mA) for weak lines. The LOX can thus be used to rank 
metal lines in terms of strength, and it can be directly compared to the detection limit in spectra 
of a given quality. Computation of the LOX as a function of [Z/H] and A'hi requires a model 
of the densities and temperatures within Lya forest absorbers, which together with the ionizing 
background radiation field determine the ionization term x'^^/x^ as a function of A^hi- 



Z 
H 



+ log(/A2) +log 



nz 



+ log 



(3) 



3. LOX for 2249 candidate metal lines 

We have applied the LOX technique to two different cosmological simulations. One is the 
standard CDM simulation described in Dave et al. (1996). A box of length 22.222 comoving Mpc 
containing baryonic and dark matter is evolved from z = 49 to 2 = 2 using TreeSPH (Hernquist & 
Katz 1989), assuming a CDM spectrum of density fluctuations with the parameters crg^-i = 0.7, 
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Q = 1.0, h = 0.5, Ofe/i^ = 0.0125, and n = 1, where h = Hq/IOO km/s/Mpc and n is the slope of the 
initial perturbation spectrum. The other simulation, a flat, nonzero-A model, also assumes a CDM 
spectrum of density fluctuations, with ag^-i = 0.8, Q = 0.4, A = 0.6, h = 0.65, Ofe/i^ = 0.0125, and 
n = 0.93. The standard model is inconsistent with the amplitude of CMB anisotropy measured by 
COBE (Smoot et al. 1992), but the normalization of the A model (computed using the CMBFAST 
program of Seljak & Zaldarriaga 1996) is matched to the 4-year COBE data. For further details 
regarding the simulation method, see Katz, Weinberg, &; Hernquist (1996a). 

The analysis of the two cosmological models yielded very similar results. The LOX is a local 
property, depending on the conditions within individual density peaks along the line of sight, 
which are insensitive to the underlying cosmology. Hence, in the following we will only describe 
the calculations for one of the models, the ACDM model, and we will focus on results at z = 3. A 
spatially uniform photoionizing background radiation field is imposed, with the shape computed 
by Haardt & Madau (1997, in preparation). They calculate the modification due to absorption 
and re-emission by the IGM of an input radiation field with J^, oc u' -i-s for A < 1050A, oc u 
for 1050A < A < 2200A, and oc u'^-^ for A > 2200A. This radiation field, at z = 3, has a larger 
drop above the He"'' absorption edge at 4 Rd than the previously published Haardt &: Madau 
spectrum (1996), by roughly a factor of three, and it is therefore in better agreement with the 
relatively large number of SilV systems observed at that redshift (see the discussion in Hellsten et 
al. 1997). 

To evaluate the ionization correction term in expression (|3|) we need to know the typical 
volume densities and temperatures in the absorbers for different values of A^hi- To find the 
relation between these quantities, we generated artificial Lya absorption spectra along 480 lines of 
sight through the simulation box at z = 3, using the method described in Hellsten et al. (1997). 
The lines of sight were selected to include a relatively large number of the high-column density 
absorbers in order to get good statistics for these rather uncommon systems. We use artificial 
spectra along random lines of sight to fix the overall amplitude of the background radiation field, 
scaling J^, by a constant factor so that the mean Lya flux decrement in the simulated spectra 
matches the value Da = 0.36 found by Press, Rybicki, & Schneider (1993; see further discussion 
by Ranch et al. 1997). 

Figure || shows results from an analysis of 1456 HI absorption features detected in these 
artificial spectra. Lines are identified by a threshold criterion in the following way. First, we define 
the velocity interval from a position where the HI optical depth increases above a critical value 
Tc = 1.2 until it drops below that value again as belonging to an HI absorption feature. Next, we 
augment this interval on both sides until either < 0.1 or a local maximum in is encountered. 
This procedure separates weakly blended "W-shaped" absorption features, if < Tc at the central 
local minimum. 

The HI column density of the absorber is found by integration of over the velocity interval, 

i.e. 
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^Hl = "^f. I r„(z;)dz; = 7.44 X 10"cm-2 / T^{v)dv, (4) 

T^e Ja-^a JAv JAv 

where v is measured in km/s in the rightmost expression. The crucial question for our purposes is 
the relation between A'^hi the typical volume density of gas in the absorbing region. Because of 
peculiar motions and Doppler broadening, the absorbing HI atoms in a given velocity bin originate 
at slightly different spatial locations. The code we developed to generate the artificial spectra 
also determines the mean total hydrogen volume density n}i{v) in the regions from which the HI 
atoms originate. In Figure |^ we plot nn vs. A'^hi, where nu is the maximum value of n}i{v) in the 
velocity interval Av of that absorption feature. In the spectra, the location of this density peak is 
close to the location of the maximum of Ta{v), as one would expect. 

As can be seen in Figure |l|, there is a well-defined correlation between the HI column density 
of a LYFA and the typical total volume hydrogen density of the absorbing region. The thick solid 
line is the power-law fit 



lognH = -14.8 + log— ^ + 0.71ogArHi, (5) 
U.U125 

which we will use to evaluate the LOX. 

Relation (^) depends on the value of Q^h"^, i.e., the total mean baryonic density. Increasing 
the value of rif,/i^ will shift the curve upward if A'hi is held constant (to satisfy the mean flux 
decrement constraint) by adjusting the intensity of the UV radiation field. The indicated scaling 
with Jlfo/i^ in (|5|) is only approximate, because it neglects effects of the temperature dependence 
on rjfc/i^ (Croft et. al. 1997). We will assume Qbh'^ = 0.0125 throughout this paper, except for the 
discussion in Section 4. While a similar correlation between nn and A'^hi holds at other redshifts, 
the constant of proportionality is different because of the evolution in the radiation field and the 
IGM structure, in particular the dilution of the cosmic mean density by the expansion of the 
universe. The nH-A^m correlation for the SCDM model is very similar, with a slope of about 0.73 
instead of 0.7, and the resulting line observability indices are nearly identical for the two relations. 

Also shown in Figure |l] is the mean total hydrogen density at z = 3 in this model. It is seen 
that the model predicts that LYFAs with log (A'^hi) ^ 14 at z = 3 arise from underdense regions 
(similar to the result of Zhang et al. 1997). 

The relation between density and temperature is slightly more complex. Figure ^ shows a plot 
of nn versus mean temperature for the same density peaks as those shown in Figure ||. For low 
densities a large fraction of the peaks follow the relation T oc n^®^. At densities nn ^ 10~^cm~^, 
the density-temperature correlation arises mainly from the interplay between photoionization 
heating and adiabiatic cooling due to the expansion of the universe. Low density regions have 
lower neutral fractions, hence lower photoionization heating rates, and they expand more rapidly, 
thus experiencing stronger adiabatic cooling. Some of the peaks lie above this temperature-density 
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relation, at T ~ lO^-^-^K, because of moderate shock heating. Peaks with lognn ^ —3.6 have 
large enough densities for radiative coohng to become significant, and the temperature is roughly 
determined by the equilibrium between radiative cooling and photoionization heating. (For a more 
detailed discussion of the density-temperature relation in cosmological simulations, see, e.g.. Croft 
et al. 1997; Hui &: Gnedin 1997; Zhang et al. 1997). For our LOX calculations we will adopt 
the relation defined by the solid line in Figure ^. This relation is again very similar to that in 
the SCDM model, although the latter shows a slightly larger scatter at lower densities. Because 
photoionization is the main source of ionization at such low densities, the LOX is insensitive to 
the exact value of the temperature. 

For an arbitrary metal line Z\ j associated with a LYFA with neutral column density A^hi 
and metallicity [Z/H], we can now evaluate the ionization correction term in (^) using the 
photoionization code CLOUDY 90 (Ferland 1996). The inputs to CLOUDY are the hydrogen 
volume density and the temperatures from the relations discussed above, as well as the UV 
radiation field. For systems with logA^ni ~ 16, continuum Hell absorption affects the radiation 
field in the central regions of the absorbers. We include this effect in an approximate way, by 
solving self-consistently for the radiative transfer through a plane parallel slab with HI column 
density 2 x A^hi, illuminated from both sides. We then use the volume averaged values of the 
ionization fractions, which should roughly be equal to those in the central region of an absorber 
with HI column density Ahi- Including this Hell-shielding effect only changes the predicted LOX 
values slightly. 

The results in Table 1 are produced by evaluating the LOX for the 2249 transitions in the 
compilation by Verner, Barthel, and Tytler (1994) and then sorting the lines according to the 
observability index. The strongest lines are shown for five different values of the HI column 
density, logA^m = 13, 14, 15, 16, and 17. We show, for logA^m > 14, only the ten strongest 
lines with A,. < 912A. These extreme ultraviolet lines are numerous, but in spectra of the faint 
QSOs at z ^ 3 they are extremely difficult to detect, even with state of the art instruments such 
as the Space Telecope Imaging Spectrograph, installed on the Hubble Space Telescope (HST) in 
February 1997 (e.g. Danks et al. 1996). At somewhat lower redshifts, the prospects should be 
better for observing these lines in STIS spectra of a few very bright QSOs, such as HE 2347+4342 
(2=2.88, Reimers et al. 1997), but the observations have yet to be made. With the Cosmic Origins 
Spectrograph (COS), scheduled to be installed on the HST in 2002 (Green et al. 1997), detection 
of these important lines may finally become possible for a larger number of QSOs at high redshift. 

A metallicity of [0/H] = —2.5 has been assumed in Table 1, and we use the solar metallicity 
values from Grevesse & Anders (1989). We evaluate the LOX for solar relative abundances, i.e. 
[Z/O]=0. It can then be scaled to any other relative abundance pattern a posteriori. For example, 
if [O/C]=0.5, as observed in some low metallicity systems (e.g. Wheeler et al. 1989), we have to 
add 0.5 to the LOX values for O lines, (i.e. use [0/H] =-2) in order to compare them to the C lines 
for [C/H]=-2.5. 
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Before discussing the results in Table 1, let us estimate detection limits in present day, high 
signal-to- noise ratio (S/N) HIRES spectra. For a Hne with rest frame wavelength A^, the rest 
equivalent width threshold for an A^^^-sigma detection is 



where A^pix = 7 is the number of pixels used for the equivalent width detection limit and Avpi^ 
(AApix) is the velocity (wavelength) spread per pixel (Churchill 1997). For the spectrum of 
Q1422-F231 obtained by SC, AA = 0.06A and Zabs « 3. Hence lines with LOX k. (2.3 - log S/N) 
should be detectable at the 5-sigma level in that spectrum. (Recall from §2 that Wr = 1mA x lO'"^"^ 
for weak lines.) The S/N per pixel is roughly 110 outside the Lya forest and 45 within the forest, 
which translates into a minimum detectable LOX of approximately 0.26 outside the forest and 
0.65 within the forest. The integration time for this spectrum was 8.2 hrs. If the QSO is observed 
for 3-4 nights the S/N ratio could be roughly doubled, so a practical absolute lower limit for 
observable lines would be LOX ^ 0. 

Looking at Table 1, the OVI(1032,1038) doublet is predicted to be a dominant metal line 
feature associated with LYFAs with logiVHi = 13, 14, and 15. For [0/H]=-2, OVI(1032) becomes 
observable for logiVni ^ 13 and tops the list for more than two orders of magnitude in Nuj. 
Because this OVI doublet has rest wavelengths less than 1216 A, the absorption features will 
always be located in the Lya forest where they are subject to heavy blending and blanketing by 
HI lines, making them difficult to identify, especially in high redshift spectra. 

A family of strong absorption lines from lower ionization states of oxygen is seen to be present 
for the higher column density LYFAs. These lines all have A,. < 833A, and cannot (for z ^ 4 — 5) 
be observed in the optical. They present a challenge for the next generation of ultraviolet 
spectrographs. 

For lines with A^ > 1216A, the CIV(1548,1551) doublet is seen to be dominant over the 
entire range of HI column densities, a remarkable result, given that this range in A^hi spans three 
orders of magnitude in volume density and hence photoionization parameter. Only for HI column 
densities close to that of Lyman limit systems (logiVHi ^ 17.2) do other metal lines, such as the 
SiIV(1394,1403) doublet, attain strengths comparable to that of CIV. This prediction is consistent 
with the fact that the first metal lines found to be associated with LYFAs were CIV lines, and 
that CIV systems are the best studied and most numerously observed metal line systems to date. 

The lines with A^ > 912A reaching the largest LOX values are CIII(977) and SiIII(1206). In 
the high column density LYFAs these lines are predicted to be so strong that the prospects of 
detecting them in the Lya forest should be reasonable, and they have indeed been observed in a 
few LYFAs in the Q1422 spectrum (SC 96). These lines deserve a more thorough investigation 
with artificial absorption spectra. A good way of testing simulation predictions for high column 
density LYFAs would be to compare CIII absorption in real spectra to the CIII features in a model 



S/N 



(Avpix/c) = 



S/N(l + Zabs) 



(6) 
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with [C/H] determined observationally from the corresponding CIV features. 

Figures 3 and 4 show graphical representations of the A'^hi-LOX relations for 12 lines selected 
from Table 1. In low density regions, the LOX is insensitive to the scatter in gas temperature, 
because photoionization is much more important than collisional ionization and because the 
recombination coefficients are not very temperature sensitive. A more complicated issue is what 
would happen if the mean temperature of the low density gas were different in the model, for 
example because of a different assumption about the heat input during the reionization period 
(Hui & Gnedin 1997; Ranch et al. 1997). The effect on the LOX is small in this case also, as we 
will argue in detail in Section 4. 

Figure 3 contains lines from three different ionization states of carbon and silicon. The 
CIV(1548), CIII(997), and CII(1335) lines are expected to be observable in absorbers with 
logiVni ^ 14, 14.9, and 15.8, respectively. These values are obtained assuming a resolution 
and S/N similar to that of SC's spectrum of Q1422+231, and a metallicity of [0/H]=-2.5. It 
is straightforward to predict these values for spectra of other resolutions and signal-to-noise 
ratios by means of equation (^, and for other metallicities by shifting the LOX-curves upwards 
by an amount A(LOX) = [0/H] -1-2.5 (or alternatively sliding the horizontal detection limit 
lines downwards by the same amount). For example, for CIV to be detectable in a LYFA with 
A'^Hi = 10^^ cm~^ in the SC spectrum, one can estimate from Figure 3 that [0/H] would have to 
be > -0.8. 

The diagonal arrows on Figure 3 indicate how the curves will shift if Ji, is multiplied by a 
factor of 2 (solid arrow) or if is changed from 0.0125 to 0.02 and Ji, is adjusted to maintain 
the observed Da (dashed arrow). The direction and magnitude of these arrows can be derived 
from approximate scaling laws, as will be discussed in Section 4. 

The modeling of the ionization correction term in the LOX is based on a single-density 
model, assuming that most of the absorption takes place at or near the maximum of the local 
density peak associated with a LYFA. This is certainly true for HI. As long as the volume density 
of a species increases with total density (i.e. the LOX increases with A'^hi) this assumption 
holds. In situations where decreases with total density this approximation becomes worse. For 
example, in our artificial spectra we sometimes see that the CII(1335) feature associated with a 
high column density absorber (logA^Hi ~ 16.5) is centered on the same velocity as HI, whereas 
the corresponding CIV (or OVI) absorption is seen as shoulders that are offset relative to the HI 
feature. This occurs because the CIV density profile does not follow that of the total density. 
Therefore, our maximum density assumption causes us to underestimate the LOX values on the 
decreasing parts of the curves in Figures 3-4. They are lower limits to the true LOX values, and 
it is necessary to do a full generation and analysis of artificial spectra if one wants to use, say, 
observed CII/CIV ratios to infer ionization properties of the IGM. Direct comparison to the LOX 
values in Figure 3 would overestimate the photoionization parameter in the central region of the 
LYFA. 
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Figure 3 also shows the sihcon lines SiIV(1394), SiIII(1206), and SiII(1260). Here, the SilV 
and Silll features are observable only for logA^ni ~ 15.8, and the Sill feature for logA^Hi ~ 16.3. 
In addition, the Sill line falls redward of the Lya forest only for LYFAs close to the Lyct emission 
peak. Clearly, for most LYFAs these silicon lines are not as promising for probing ionization 
conditions as the carbon lines discussed above. 

Lines from the family of very strong O transitions 0111(833), OIV(788), OV(630), and 
OVI(1032), as well as NV(1239) and SVI(1063), are displayed in Figure 4. OVI is by far 
the strongest 'high ionization' line with > 912A, and it is observable in the entire interval 
13.6 ^ logATni ^ 15.8. For systems with logTVHi ~ 14.5, OVI should be detectable in the SC 
spectrum, provided [0/H] ^ —3. In our model, the OVI is produced mainly from photoionization, 
and the LOX is not very sensitive to the assumed temperature. With present day spectra, OVI is 
the best and only candidate for probing the metallicity of the IGM in regions close to the mean 
density. As we already mentioned, observations of the OVI lines are hampered by the high line 
density within the Lya forest. Low column density lines in cosmological simulations are usually 
broadened by Hubble flow and peculiar velocities rather than thermal motions, so the OVI lines 
may not necessarily be distinguishable by low 6-parameters. In individual cases, however, it should 
at least be possible to place upper limits to the metallicity of the systems that happen not to have 
strong HI contamination at the expected positions of the OVI features. 

The NV(1239,1243) is predicted to be weaker than the OVI doublet and probe a somewhat 
denser part of the IGM. If [N/C] is as low as -0.7 as suggested by Pettini et al. (1995) the NV 
doublet will be observable only in LYFAs with relatively high metallicity, roughly [0/H] ^ —1.8. 
This seems consistent with the fact that SC observed NV in only one out of six candidate absorbers 
in their Ly/3 selected sample. 

Towards lower redshifts, where the Lya forest begins to thin out, it should become increasingly 
easy to identify OVI absorption lines by their coincidence with corresponding Lya and Ly/3 
absorption, at least in a statistical way. Though somewhat arduous, this type of investigation with 
OVI is probably the most promising route to testing the trend between metallicity and HI column 
density predicted by some models of early IGM enrichment (e.g., Gnedin &: Ostriker 1997). 

4. LOX dependence on model parameters 

Figures 3 and 4 present LOX results for a single redshift and a single set of model parameters: 
ACDM at z = 3, with fi^/i^ = 0.0125 and the Haardt-Madau (1997) UV background spectrum 
scaled in intensity to give an HI photoionization rate Fhi such that Da = 0.36 in the simulated 
spectra. However, we can estimate the effects of changing some of these parameters by simple 
scaling arguments because the dynamics of the gas in the density regime relevant to LYFAs is 
driven by the gravity of the dark matter, making the spatial distribution of the gas insensitive to 
changes in O;,, Jv, or the IGM temperature. We can therefore compute the effects of such changes 
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on the A'^Hi vs. riH correlation (Figure |l|), and we can compute the effect on LOX values through 
equation (^) if we make the approximation that the ionization fractions depend only on the 
ionization parameter T/rin. Because recombination rates vary only slowly with temperature and 
the range of temperatures in LYFAs is much smaller than the range of densities, this is a good, 
but not a perfect, approximation. 

Consider first the effect of changing Ju by a constant factor q. This change in isolation would 
alter the mean flux decrement Da, which we have thus far adopted as an observational constraint 
on the model, but current observational uncertainties in allow some range in Ju (compare, 
e.g., Zuo & Lu 1993 to Press et al. 1993 and Ranch et al. 1997), and comparison of simulated 
and observed HI column density distributions favors a somewhat higher value of Ju than we 
have adopted here (Dave et al. 1997; Gnedin 1997a). We assume that the shape of J^ does not 
change, so that the photoionization rates T for all species change by the same factor. An absorber 
that previously had column density iV^j now has a column density A^'hi = g~^iV^j — raising J^ 
lowers the hydrogen neutral fraction and therefore lowers A'^hi- The density nu in this absorber 
is the same as before, so the ionization parameters T/nn are the same as those in an absorber 
that previously had column density A^^j = Q'~^'^^A^jji ~ 'i'~'^'^^-^Hi- We have used the correlation 
nn oc A^Hi'^'^ from Figure |l| to convert a factor change in nn to a factor q~^-^^ change in A^hi- 
Assuming that log(x^/2;^) depends only on T/nu, we see from equation (^) that 

LOX(gJ^,logA'Hi) = LOX(J^, log AThi -0.43 log g) + 0.43 log g, (7) 

where the second term accounts for the difference logA^Hi — log A^^j. For example, the effect of 
doubling J^ while keeping other parameters fixed is to shift the LOX(log A^hi) curves diagonally 
by the vector (0.13,0.13), indicated by the solid arrow in Figure 3. The magnitude of this shift is 
small, and the direction is nearly parallel to the LOX curve itself. The same shift would apply to 
Figure 4. 

If we keep J^ fixed but change 0^ by a factor q, then absorber densities tt-h change by 
q and column densities by q^ (because the hydrogen neutral fraction is itself proportional 
to tt-h). An absorber that previously had column density A^^j now has A^Hi = q^%i- The 
ionization parameters are the same as those in an absorber that previously had column density 
^Hi ~ Q^'^^^m — 9 "■^''^Hi- With the same reasoning as before, we find 

LOXiqQb, log A'hi) = LOXi^b, log A^m + 0.57 log q) - 0.57 log q, fixed J^. (8) 

Since the observational constraints on Da are better than those on the ionizing background 
intensity, a factor q change in 0;, should really be accompanied by a factor change in J^, so 
that Da and the HI column densities of absorbers remain the same. The net change is to increase 
r/nn by a factor of q at fixed A^hi- An absorber with column density A^hi has the same ionization 
parameters as an absorber that previously had column density A^^j = (7~^'^^A^hi- Thus, 



LOXiqnb,q^Ju,logNm) = LOX(Qfe, J^, log A-ri - 1.43 log g) + 1.43 log g, fixed Da- (9) 
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For example, the dashed arrow in Figure 3 shows the effect of increasing Qb from the value 
Qb = 0.03 used in our simulation to the value Qb = 0.047 favored (for h = 0.65) by the low [D/H] 
value of Tytler et al. (1996). 

Plausible variations in the assumed gas reionization history can alter the temperature of the 
low density IGM by factors ~ 2 — 3 (Miralda-Escude & Rees 1994; Hui & Gnedin 1997). Because 
the HI recombination coefficients are cx T^^-"^ for T ~ 10^ K, a factor q change in gas temperatures 
requires a factor change in Jy to keep Da fixed. From the scaling relation (0), we see that 
the impact of such a temperature change would be small. 

To estimate the evolution of the LOX with redshift, we can appeal to the argument of 
Hernquist et al. (1996) and Miralda-Escude et al. (1996) that the primary driver of evolution in 
the Lya forest at high redshift is the expansion of the universe, which changes the cosmic mean 
density in proportion to (1 + z)"^ . Non- linear gravitational evolution changes the overdensity 
field p/p, but this evolution has a smaller impact on the forest than the overall expansion, at 
least for low column density absorbers. Matching the observed evolution of Da requires an HI 
photoionization rate Fhi that is roughly constant between z = 4 and z = 2 (Hernquist et al. 
1996; Miralda-Escude et al. 1996; Ranch et al. 1997), so if we ignore changes to the shape of the 
spectrum we can take Jy to be constant over this interval. A factor of two decrease in (1 + z)^ — 
roughly the change between z = A and z = 3 or between z = 3 and z = 2 — has the same effect as 
a factor of two ZTlCTCdSC 111 Ji/ clt a given redshift, illustrated by the solid arrow in Figure 3. While 
this argument is only approximate, it is enough to show that the LOX should not vary strongly 
with redshift over the range z = 2 — 4, except in the case of ionic species that are affected by the 
change in the spectral shape of Jjy. 

Structure formation models that have a power spectrum amplitude similar to that of ACDM 
on Mpc scales at z = 3 will predict similar spatial and density structure in the diffuse IGM. 
(Models with very different power spectra will of course predict different structure, but such 
models may soon be ruled out by the Lya forest data themselves [see, e.g., discussions by Ranch 
et al. 1997; Gnedin 1997b; Croft et al. 1998].) If the baryon density ^bh"^ is held fixed, then the 
most important effect of changing the cosmological model is to alter the expansion rate H{z). 
Matching Da requires Fhi oc H~^{z) (see, e.g., equation [2] of Ranch et al. 1997), so halving H{z) 
has the same impact as doubling J^. In practice the difference between H(z) values in currently 
popular models at z ~ 3 is less than a factor of two, so the expected impact is small. 

We conclude that our quantitative predictions for the LOX are likely to hold rather generally 
in the cosmological picture of LYFAs, at least for column densities A'^hi ^ 10^^ cm^^, where 
the scaling arguments used in this Section are most secure. Our comparison of two numerical 
simulation models, a standard CDM model and a ACDM model, supports this claim, and it does 
not show any major changes to the LOX predictions even for higher column densities. Of course, 
the LOX is always directly affected by the metallicity of the absorbing gas (equation Q), and 
the insensitivity to cosmology implies that absorption observations can therefore yield robust 
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estimates of the metallicity in different regions of the IGM. Ionization corrections are essential in 
the interpretation of such measurements, but the corrections are not strongly dependent on the 
assumed cosmological model. 

5. Discussion 

In this paper we have introduced the LOX, a measure that allows a comprehensive selection 
of the metal lines that are useful for testing the detailed predictions of models for LYFAs. We 
have applied this method to a particular model, a cosmological ACDM simulation, and we have 
quantitatively discussed the dependence of the results on the model parameters J,^, Clj,, and z. We 
have also performed this analysis for a different cosmological simulation, a standard CDM model, 
but the results differ negligibly from those of the ACDM model. 

We find that very few detectable metal lines are associated with the low to moderate density 
regions of the IGM that give rise to absorbers with logA^ni ~ 15- Most of the baryons in the 
Universe at z = 3 are believed to reside in this part of the IGM (e.g., Zhang et al. 1997), and the 
OVI(1032,1038) doublet is found to be the best probe of the metallicity in these regions. Other 
strong O lines, such as OV(630), OIV(788), and 0111(833) are technically difficult to observe, and 
have yet to be detected. The CIII(977) and SiIII(1206) lines are the strongest metal lines in high 
column density LYFAs. The CIV(1548) line is the strongest metal line with > 1216A for all 
LYFAs with log A^hi ^ 17. While the specific values of the LOX as a function of column density 
will depend on redshift and on the adopted cosmological model and radiation field, we have argued 
that the qualitative conclusions presented above are insensitive to the choice of cosmological model 
of high redshift LYFAs, and they should hold generally within this scenario. 

The relatively uncommon absorbers with logA^Hi ^ 16 should offer the strongest test of 
cosmological simulations. These systems should have more than a dozen observable metal lines 
associated with them (for [0/H] ^ —2.5). The relative strengths of these lines will depend 
on the assumed radiation field and abundance pattern, so a first check would be whether it is 
possible to make the line strengths in the model spectra match those of lines associated with 
high column density LYFAs for reasonable values of model parameters. Then one should make a 
detailed examination of the relative locations in velocity space between absorption components 
of individual species. A detailed comparison to observations will ultimately determine if the 
resolution of the simulations is adequate in these high-density regions, and to what extent local 
dynamical effects associated with the production of metals can be neglected. There are some hints 
of discrepancies between models and observations for these high column density systems. For 
example, the models in some cases predict the velocity components of CIV and CII absorption 
lines to differ, whereas there seems to be some new observational evidence that CII and CIV 
components are found at the same positions (A. Boksenberg, private communication). 

Cosmological simulations depict the typical low column density LYFAs as simple structures: 
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low density, smooth, and governed by the straightforward physics of gravity, cosmic expansion, 
and photoionization. Studies of absorption along parallel lines of sight provide direct observational 
support for this point of view, especially recent HIRES observations of gravitationally Icnsed QSOs 
that set stringent lower limits to the scale of any substructure within low column density absorbers 
(Ranch 1997; see also Smette et al. 1992, 1995). The simulations do not predict a sudden change 
in LYFA properties at any particular column density, but they do predict a steady increase of 
gas density with A?hi (Fig. 1). For logiVni ^ 16, the typical absorbers are partially collapsed 
structures, and gas within them has begun to experience radiative cooling. In this regime, the 
finite mass and spatial resolution of the simulations may limit the accuracy of the results, and 
local astrophysical processes such as star formation and supernova explosions might influence the 
structure of these systems. The standard CDM model studied here produces only about 1/3 of the 
observed number of Lyman limit systems (logA^Hi > 17.2), even though it matches the abundance 
of much stronger, damped Lyct absorbers quite well (Katz et al. 1996b; Gardner et al. 1997). 
While further work is needed to assess the sensitivity of this result to cosmological parameters 
— especially the baryon density parameter fib — the discrepancy of numbers suggests that a 
population of absorbers unresolved by the simulations may become important at column densities 
logATni ~ 17. 

The high column density LYFAs in the cosmological simulations represent the transition 
between the weak IGM fluctuations traced by low column density lines and the cold gas 
concentrations in virialized dark halos traced by damped Lya absorption. The high-LOX lines 
listed in Table 1 can reveal many details of the structure and physical conditions in these 
absorbers, testing the accuracy of the simulations in this regime, providing clues to the nature of 
any additional absorber populations, and capturing snapshots of gas as it makes its way from the 
IGM into high redshift galaxies. 
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Fig. 1. — Maximum total hydrogen number density in units of log (cm~^) versus the HI column 

density in units of log (cm~^) for 1456 HI absorption features along 480 lines of sight at z = 3 in 
our cosmological model. The solid line shows the fit lognn = —14.8 + O-TlogiVni. Also shown is 
the mean hydrogen density at z = 3 (dashed line). 

Fig. 2. — Temperature vs. hydrogen number density for the same density peaks as plotted in figure 
1. The shown fit is described by T = 7.17 + 0.65 lognn for logT < -3.8 and T = 3.18 - 0.41ognH 
for logT > -3.8. 

Fig. 3.— L0X(7Vhi) for CIV(1548), CIII(977), and CII(1335) (solid curves) and SiIV(1394), 
SiIII(1207), and SiII(1260) (dashed curves) at z = 3, assuming [0/H] = -2.5. The horizontal 
lines are approximate detection limits within the Lya forest (dashed line) and redwards of the 
forest (solid line) in Songaila and Cowie's spectrum of Q1422+231 (equation (6)). Arrows indicate 
the effects of doubling the radiation intensity (solid arrow) and of changing from 0.0125 to 

0.02, as discussed in Section 4. 



Fig. 4.— LOX(ArHi)for OVI(1032), OV(630), OIV(788), and 0111(833) (solid curves), and 
NV(1243), and SVI(933) (dashed curves). 
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Table 1: LOX sorted lines for five different values of log A^hi 
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